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ABSTRACT 

After a brief overview of the so-called silent models and their present status, we 
consider the subclass of Bianchi Type-I models with a magnetic field source. Due 
to the presence of the magnetic field, the initial singularity shows "oscillatory" 
features reminiscent of the Bianchi Type-IX case. The Bianchi Type-I models 
with a magnetic field are therefore a counterexample to the folklore that matter 
fields can be neglected in the vicinity of the singularity. 
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1. Introduction 



Often a preferred timelike vector field u a exists in a cosmological context (usu- 
ally the 4- velocity of the fluid), and various kinematical and dynamical quantities 
are defined with respect to u a : for example, the electric and magnetic parts of the 
Weyl tensor C a bc d are, respectively, E a b := C a c bd,u c u d and H a b := \r] a c J C e f bd u c u d . 
The so-called "silent universes" are cosmological models characterized by two spe- 
cific features: the first regards their matter content, which is taken to be a perfect 
fluid with vanishing pressure, p = 0; the second regards the dynamics, constrained 
by assuming the vanishing of the magnetic part of the Weyl tensor, H ab = 0. The 
name "silent" itself comes from noticing that the assumptions p = H ab = imply the 
lack of sound and gravitational wavesSJ In the recent cosmological literature, models 
characterized by p = H ab = plus a third assumption, the vanishing of the vorticity 
of the fluid oj ab = 0, were originally introduced by Matarrese et al (1993) to model the 
non-linear evolution of perturbations during the matter dominated era.B It was hoped 
that these models could give a reasonably good description of the growth of struc- 
ture in a collisionless medium. The assumptions p = uo ab = are very standard in 
dealing with these problems; the restriction H ab = was introduced by extrapolating 
from the vanishing of EL b in first-order scalar perturbations.B Matarrese et al (1993) 
also implicitly assumedcl that the field equations were compatible with the H ab = 
constraint or, in other words, that the evolution of initial data preserves H ab = 0. 
This assumption was corroborated by a preceding paper by Barnes and Rowlingsoni 
(1989), and later on by work! of Lesame et al (1995). 

In studying models with p = u ab = H ab = one takes the matter density /i, the ex- 
pansion scalar 0, the shear a ab and the electric Weyl tensor E ab as physical variables; 
the dynamics of these quantities is then determined by a set of field equations arising 
from the Ricci and the Bianchi identities? The amazing feature of these models is that 
the evolution equations for the variables /i, 6, a ab , E ab become a set of quasi-linear 
first-order ordinary differential equations (ODEs), as was first explicitly showncl by 
Matarrese et al (1993). This fact itself was enough to raise the attention of cosmolo- 
gists on these silent modelsH would they provide a fair description of the non-linear 
evolution of cosmological perturbations? The assumptions p = uj ab = H ab = would 
simplify computer simulations of large-scale structure formation enormously. Unfor- 
tunately, it was soon realized by Matarrese et ala (1994) that the H ab = assumption 
already breaks down at the second order in the perturbations. In this context it should 
be noted that the cosmological effects of H ab have never been properly quantified. As- 
suming that the constraint H ab = was compatible with the evolution equations for 
the physical variables, the dynamics of silent models with p = uo ab = H ab = was 
investigated from a dynamical system point of view by Bruni et al (1995), who fo- 
cused on the character of the singularities.0 They found that the generic final fate of 
a fluid element under the p = u ab = H ab = restrictions is to fall into a spindle-like 

*In this approach, the Einstein equations are used to establish an algebraic constraint between the 
Ricci part of the curvature, and the energy-momentum tensor. 



Kasner singularity, while only a set of measure zero of initial data leads to a pancake 
(which is instead thought to be generic in a Newtonian collapse scenario). Finally, 
in a parallel work, Bruni et al (1995) included a cosmological constant A in their 
analysis^ in order to investigate the so-called "cosmic no-hair" conjecture,^ and it 
was indeed found that the de Sitter spacetime is the unique attractor for initial data 
leading to expansion! An important class of spatially inhomogeneous exact solutions 
resides in the silent family, namely the well-known Szekeres modelsHj These contain 
the spherically symmetric Lemaitre-Tolman-Bondi models which in turn have the 
spatially homogeneous Kantowski-Sachs models as a special subclass. The Szekeres 
models are of Petrov type D, and may be called degenerate, in that the Weyl tensor 
(reducing to its electric part E ab ) and the shear a ab each have two equal eigenvalues. 
Recently the existence of a more general spatially inhomogeneous class of silent mod- 
els (of non-degenerate Petrov type I) has been challenged.0lll Furthermore, even if 
these models exist, their singularities are unstable agains t H* perturbation^ and 
their relevance to structure formation is at least doubtful.El Therefore, in the following 
we shall focus on homogeneous spatially flat silent models of Bianchi Type-I with a 
magnetic field, which do exist as solutions of the Einstein-Maxwell field equations. 

2. Dynamics of spatially flat silent universes with a magnetic field 

As we said above, the evolution equations one obtains under the assumptions 
p = uj a b = H a b = are ODEs. Then, it is natural to ask if these restrictions can 
be relaxed, or if other matter fields can be considered, such that the ODE-structure 
of the evolution equations is retained. We have found that this is indeed the case if 
one has non-vanishing fluid vorticity u a b ^ 0, or if one includes a magnetic field H a 
as an additional matter source. In the latter case also the Maxwell equations come 
into play. When p = u> a b = Hab = H a = 0, the shear a a b and the electric Weyl tensor 
E a b admit a common Fermi- Walker transported eigentetrad.HI This is also true if a 
magnetic field is aligned with one of the shear eigendirections. However, for non-zero 
fluid vorticity u a b or a general magnetic field H a , the eigentetrads of the shear and the 
electric Weyl tensor are no longer Fermi- Walker transported, nor do they coincide. 

One can show that non-shear-aligned magnetic fields are only compatible with a 
spatially homogeneous Bianchi Type-I spacetime geometrypa which we will consider 
in the following. The kinematical quantities we use refer to the normal timelike 
congruence associated with the spacelike surfaces of homogeneity. For simplicity we 
take the energy-momentum tensor to be comprised of the magnetic field H a only."'" The 
Bianchi Type-I spacetime geometry automatically implies uj a b = H a b = 0. Here, we 

tit was also found! that there are also initial data for which, even after a period of inflation, a 
recollapse occurs, in which case A is soon negligible, and the collapse leads to a Kasner spindle-like 
singularity, as in the A = case. On the other hand, when A = 0, the attractor for initial data 
leading to expansion is the Milne modellll 

■'■A dust matter field with pn = and/or a cosmological A can always be added, and indeed their 
coming into play makes the discussion of the corresponding models more interesting for physical 
purposes, but for the sake of simplicity we shall not consider them here. 



adopt a tetrad description, where we denote the spatial components of tensor fields 
by Greek indices, and we choose a shear eigentetrad. From a dynamical point of view, 
the configuration outlined is completely characterized by two of the shear eigenvalues, 
e.g., a 2 and 03, the expansion scalar 9, and the three magnetic field component H a . 
However, in order to study their dynamical properties, it is very useful to introduce 
five dimensionless variables, a dimensionless angular velocity (Fermi- rotation), and a 
new time derivative: 

£±:=<7±/e, H a :=V3H a /Q, R a :=n Q /Q, d :=3e /&, (1) 

where a + := 3(<7 2 + o^)j2 and cr_ := V3(a 2 — cr 3 )/2, Q a gives the angular velocity 
of the spatial tetrad with respect to a Fermi-transported basis, and eo denotes the 
time derivative along u a . Finally, it is also useful to define the standard deceleration 
parameter q : = — 39/9 2 — 1 and a density parameter Q H := H. a H a /2, related to the 
shear variables £± by 

g=l + £+ 2 + £_ 2 ; ^ = l-£ + 2 -£_ 2 . (2) 

With these definitions, we obtain the following evolution equations: 



d x+ = 


-(2 


-q) Z + + H 1 2 -±(H 2 2 + H 3 2 ) , 


(3) 


a s_ = 


"(2 


-g)S_-f (ft 2 2 -H 3 2 ), 


(4) 


doH! = 


"(1 


-g + 2S + )^ 1 + 3(i? 2 H 3 -R 3 H 2 ) , 


(5) 


d H 2 = 


"(1 


-9-E+- v / 3S_)H 2 + 3(i? 3 ^i-^i^3) , 


(6) 


d H 3 = 


-(1 


+ v / 3S-)H 3 + 3(i?iH 2 -i? 2 7i:i) , 


(7) 



where the R a are determined by the algebraic constraints 

^ 1 ^ 2 +3(S 2 -S 1 )i? 3 = 0, H 2 H3+3(S 3 -S 2 ) J R 1 = 0, H 3 ^i+3(S 1 -S 3 )i? 2 = 0, 

(8) 

where Si = — 2S + /3, S 2 = (S + + v / 3S_)/3, S 3 = (E + — y / 3S_)/3. The physical 
phase space of this dynamical system admits various invariant submanifolds: the most 
important for the present discussion are the three that arise when the magnetic field 
is aligned with one of the shear eigendirections. The latter can be Fermi-transported 
in each of these cases, e.g., 7i 2 = 7i 3 = 0, as it follows from (§) in the non degenerate 
case. Using the algebraic constraints @) one sees that now equations are 
reduced to a planar system: 

d £+ = (2-S + )(l-S + 2 -S_ 2 ), (9) 

O £_ = -S_(l-E + 2 -S_ 2 ), (10) 

where from (0) one has 1 — S + 2 — XL 2 > 0. The system above is the heart of the 
dynamics of the more general case when the magnetic field is not aligned; it has one 



isolated stationary point at S + = 2, £_ = (outside the physical phase space), and 
is also stationary on the whole Kasner ring S + 2 + XL 2 = 1. The eigenvalues of the 
Jacobian at points on the Kasner ring are Ai = and A2 = — 4£ + + 2. The first 
vanishes because there is no tangential motion along the ring, the second is non- 
negative, A2 > 0, for S + < 1/2. Thus, one third of the Kasner ring is asymptotically 
stable, because of the aligned magnetic field. These local properties of the Kasner 
ring are better understood noticing that the planar system above admits the explicit 
solutions 

£_=C(£+-2), - i/Vs < c < i/Vs , (11) 

which are heteroclinic orbits joining pairs of points on the ring (they are straight lines 
passing through the point S + = 2, E_ = 0). Since the magnetic field can be aligned 
with any of the three spatial directions, there are two other stationary points of the 
system (|3])-([7|) in the (£+, IL)-plane, located at S + = — 1, £_ = ± \/3. Passing 
through these two points, two other sets of heteroclinic orbits arise. The overall 
effect of this structure in the physical phase space is that the Kasner ring becomes 
completely unstable, as each of its points becomes a saddle. 

3. Conclusions 

The phase space structure in this example is reminiscent of the Bianchi Type- 
IX dynamics, in that the Kasner ring no longer is a source at any point. Thus, 
the approach towards the initial singularity becomes more complicated and probably 
leads to an "oscillatory" regime. Unfortunately, it is harder to obtain a detailed 
understanding of this initial behaviour in the present situation because of the existence 
of the constraints (H). We do not want to address the question if this behaviour is 
chaotic or not. This theme has been debated for a long time in the context of the 
Bianchi Type-IX models (see the contribution of Di Bari et al in this volume). Clearly, 
further work is needed to clarify the dynamics in the vicinity of the singularity of the 
Bianchi Type-I models with a magnetic field. In particular, a numerical investigation 
would be desirable. 

The Kasner ring is a source for Bianchi Type-I dust spacetimes with no magnetic 
We have seen that this is no longer the case when magnetic fields are included. 
Indeed, the important point is that a magnetic field completely changes the structure 
of the phase space in the vicinity of the initial singularity, and as one gets closer 
to the singularity this effect becomes more and more relevant. Hence, contrary to 
some folklore, matter fields can strongly influence the dynamics close to the initial 
singularity. 
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